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Mirror dark matter, and other similar dissipative dark matter candidates, need an 
energy source to stabihze dark matter halos in spiral galaxies. It has been suggested 
previously that ordinary supernovae can potentially supply the required energy. 
By matching the energy supplied to the halo from supernovae to that lost due to 
radiative cooling, we here derive an approximate scaling relation, Rsn oc PqTq {Rsn 
is the supernova rate and po, tq the dark matter central density and core radius). 
Such a relation is consistent with dark matter properties inferred from studies of 
spiral galaxies with halo masses larger than 3 x IO^^Mq. We also present a line of 
reasoning, related to the details of supernova energy transport in the halo, which 
explains why galaxies have a cored distribution and yields a second scaling relation, 
tq oc 1/po- This relation is also in agreement with observations. 
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A hidden sector exactly isomorphic to the ordinary matter sector is required if one 
hypothesizes fundamental improper space-time symmetries which are unbroken[l]. 
In such a theory, there is a mirror particle corresponding to every type of ordinary 
particle, except perhaps the graviton. Thus, a spectrum of stable dark particles 
naturally arises. We denote these mirror particles with a prime ('), e',H\He\ .... 
The symmetry implies that the masses of these particles are identical to their ordi- 
nary matter counterparts and that the mirror particles interact with mirror gauge 
fields (such as the mirror photon) in a manner completely analogous to the ordinary 
matter sector. 

Mirror particles have emerged as an interesting candidate for dark matter (for re- 
views and more complete bibliography see e.g. [2]). Mirror dark matter can explain[3] 
the positive dark matter signals from the DAMA[4], CoGeNT[5] and CRESST-II[6] 
direct detection experiments. This requires photon-mirror photon kinetic mixing 
(defined below) of strength e ~ 10~^. Mirror dark matter can also explain[7] the 
large scale structure of the Universe (matter power spectrum and CMB) in a man- 
ner analogous to standard coUisionless cold dark matter models provided [8] that 
e ~ (1 - 3) X lO-*^. 

On small scales mirror dark matter has a number of distinctive features due 
to self-interactions and dissipative interactions. In contrast to coUisionless parti- 
cles, galactic halos of spiral galaxies are composed predominately of mirror particles 
in a pressure supported spherical plasma [9]. There may also be a subcomponent 
consisting of compact objects such as old mirror stars[10] (gravitational lensing ob- 
servations limit the MACHO halo fraction of the Milky Way to be less than around 
0.3[11]). Because mirror dark matter is dissipative, an energy source is needed to 
stabilize the dark matter halos in galaxies. [Without an energy source the mirror 
particles would collapse to a dark disk on a time scale typically around a few hun- 
dred million years.] It has been speculated previously[9] that ordinary supernovae 
can potentially supply the required energy if photon - mirror photon kinetic mixing 
exists [12]: 

where F^j^ {F'^^^) is the field strength tensor for the photon (mirror photon). The 
physical effect of the kinetic mixing interaction[13] is to induce a tiny ordinary elec- 
tric charge (oc e) for the mirror charged particles. In the hot and dense core of type 
II supernovae mirror electrons and positrons can be produced from plasmon decay 
processes[14]. Thus ordinary supernovae can be a source of light mirror particles as 
well as the ordinary neutrinos. Indeed, it is estimated that half of the core collapse 
supernova energy is emitted by mirror particles (e', e', 7') if e ~ 10^^ [14, 15]. A 
significant fraction of this energy might possibly be absorbed by the mirror parti- 
cle halo. If this fraction is independent of galaxy size then an approximate scaling 
relation, Rsn oc p^r^ {Rsn is the supernova rate and po, '"o the dark matter cen- 
tral density and core radius) follows by matching the energy supplied to the halo 
from supernovae to that lost due to radiative cooling. We find that this derived 
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relation is roughly consistent with dark matter properties inferred from studies of 
spiral galaxies. We also present a line of reasoning, related to the details of super- 
nova energy transport in the halo, which potentially explains why galaxies, in this 
picture, have a cored distribution and yields a second scaling relation, tq oc 1/po- 
This relation is also in agreement with observations. Although our discussion is in 
the context of mirror dark matter, similar conclusions could be obtained for a class 
of dissipative dark matter candidates which closely mimic mirror dark matter, such 
as those discussed recently in ref.[16]. 

The physical picture then, is that galactic halos in spiral galaxies arc composed 
predominately of mirror particles e',H',He', .... in a self-interacting pressure sup- 
ported halo. Clue's about the chemical composition of this halo arise from early 
Universe cosmology. Calculations indicate[17] a primordial mirror helium mass frac- 
tion around 0.9 for e ~ 10~^. This suggests a halo composed primarily of He' with 
perhaps a fraction of H' and mirror metal components (produced in mirror star 
formation at an earlier epoch). To a first approximation, we can consider the halo 
as composed of mirror helium, which for temperatures above around 40 eV should 
be fully ionized. The radiative cooling rate, Tcooi, of such a mirror particle plasma 
is given by the analogous expression for ordinary matter plasma[18] 



where is the e' number density, and A(T') is the cooling function and has units 
of erg cnr^ s^^. For a mirror helium plasma, ng/ ~ 2p(im/^He- 

Rotation curves in spiral galaxies arc well fit [19] with a baryonic component 
described by a Freeman disk[20] and a spherically distributed cored dark matter 
component with the Burkert profile [21]: 



where tq, po is the dark matter core radius and central density respectively. However, 
as discussed in ref.[9], a spherical self gravitating isothermal gas of particles requires 
p oc 1/r^. If the temperature is not isothermal but rises in the inner region then 
it is possible to produce a cored distribution. This motivates the quasi-isothermal 
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profile with p^m — i '^l^-i^ ■ Note that if, instead of the Burkert profile we adopted the 

quasi-isothermal profile, which also provides a reasonable fit to the rotation curves of 
spiral galaxies, this would not significantly affect our subsequent analysis. Anyway, 
assuming a mirror dark matter halo composed (predominately) of an ionized plasma 
with the Burkert density profile, it follows that 




(2) 



Pdm 



(3) 





(4) 



where 
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Thus, we find: 

- (10-23 ergim^/s) [w^t^J^] (lOl^) ' 

On the other hand, the rate at which energy is supphed to the halo by ordinary core 
collapse supernovae is: 

^SN — fsN{EsN)RsN 

where fsN is the proportion of the supernova total energy, Esn, supplied to the halo 
and RsN is the galactic supernova rate. Equating Tcooi with Vsn for the Milky Way 
galaxy implies fsN ~ 0.1 with an order of magnitude uncertainty. 

In this picture, Tcooi — '^sn should hold for any galaxy, not just the Milky Way. 
Imposing this condition suggests a scaling relation: 

RsN oc A(r)pgrg . (8) 

The idea is that the dynamics can keep this relation satisfied. If Tcooi > ^sn {^cooi < 
Tsn) then the halo should contract (expand), thereby increasing (decreasing) the 
star formation rate, and hence Rsn, until Tcooi ~ Tsat. What is the temperature T? 
For an isothermal halo in hydrostatic equilibrium, we expect [9] 

T ^ \mvl^, (9) 

where fa is the mean mass of the particles in the halo (m i=s 1 GeV for a mirror 
helium dominated halo) and v^ot is the asymptotic value of the rotational velocity. 
For spiral galaxies, halo masses have values 3 x 10^° Mq ~ Mh ~ 3 x 10^^ Mq. 
For such halo masses, v^ot has the range 50 km/s ~ Vrot ~ 500 km/s, and Eq.(9) 
suggests a rough temperature range: 

10 eV ~ T ~ 1000 eV . (10) 

At the lowest temperature's mirror helium will not be fully ionized and the cooling 
rate can become dramatically suppressed. This might explain why low mass halos 
< 10^° Mq hosting disk systems are not detected. For temperatures above ~ 40 eV, 
mirror hehum is fully ionized and for 40 eV ~ T ~ 1000 eV, A(T) typically varies by 
only a factor of 2-3 [18]. [This is relatively minor compared to the variation of, say, 
ro over this range of halo masses, which is around 2 orders of magnitude.] Taking 
the approximation of A(T) as constant in this temperature range, yields the scaling 
relation: 

Rsn oc plrl . (11) 
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Is the above relation satisfied by spiral galaxies? A set of scaling relations have 
been derived from observations of rotation curves in spiral galaxics[22, 23, 24, 25] 
and summarized recently in ref.[26]. These relate Po,rQ, and r*-band luminosity, 

~ 2.2 ±0.25 




:5xiOii,U 



2.65 




1 + 



3xlO"M0 



2.00 



~ 0.66 + 0.58 log ( | . (12) 

These relations imply that for spiral galaxies with ~ 3 x IO^^Mq, 

Lr oc rl'^ . (13) 

Observations [27] indicate that the galactic B-band luminosity. Lb-, scales wtih the 
type II supernova rate as oc R^s%i with an uncertainty in the exponent around 
0.1. Neglecting the expected minor difference between the scaling of and we 
arrive at an approximate 'empirical' scaling relation: 

RsN « rS-« . (14) 

This relation is consistent, within the uncertainties, to the rough 'theoretical' scal- 
ing relation arrived at in Eq.(ll) (given that po^o is observed to be approximately 
constant). This provides some observational support to the notion that ordinary 
supernova supply the energy needed to stabilize halos in spiral galaxies, at least for 
Mfe ~ 3 X W^Mq. What about for ~ 3 x lO^^M©? These halos are expected 
to have a lower temperature T 40 cV. for which the free mirror electron fraction 
reduces as e' become bound into atoms. This means that the right-hand side of 
Eq.(ll) will begin to reduce more rapidly as a function of tq. This feature appears 
quahtatively consistent with observations as indicated by the sharper fall of the 
right-hand side of the second equation in Eq.(12) for Mh ~ 3 x lO^^M©. 

The above analysis has assumed that a significant fraction {fsN ~ 0.03 or so) of 
supernova energy can be transmitted to the halo and that this fraction is roughly 
independent of galaxy size. How reasonable arc these assumptions? Let us assume 
a kinetic mixing parameter e ~ 10~^, so that around half of type II supernova 
energy is converted into e', e! emitted from the core with energies ~ MeV. One 
could imagine that the huge number of MeV e', e' injected into a volume [(~ 1 pc)^] 
around ordinary supernova will radiatively cool, converting most of their energy 
into mirror photons. The energy spectrum of these mirror photons is of course very 
hard to predict but it might be roughly analogous to the 7 spectrum of ordinary 
Gamma Ray Bursts (GRB's). GRB's feature a wide spectrum of energies with mean 
around 700 keV with a few percent of energy radiated below 10 keV. In any case. 
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these mirror photons will then heat the mirror particle halo, potentially supplying 
the energy lost from the halo due to radiative cooling. Whether this can happen 
depends on how strongly the mirror photons scatter off mirror electrons, both bound 
and free. 

Consider first the scattering off free mirror electrons, i.e. elastic (Thomson) 
scattering with E'^ independent cross-section ar — 6.7 x 10~^^ cm^. We estimate 
that the optical depth due to elastic scattering for 7' propagating out from the 
galactic center is 



POO 

tes = / arUe'dr 
Jo 



^ 0.78crrporo ( ^] ~ 0.006 (15) 

where wc obtained PqTq from Eq.(12). We expect /s7v~t and Eq.(7) then suggests 
that elastic scattering is probably not frequent enough to supply enough heat to the 
halo to stabihze it. The cross-section is around an order of magnitude too small. 
However if the halo contains a significant proportion of heavy mirror elements - 
necessary to explain the direct detection experiments [3] - then the photoelectric 
cross-section of heavy mirror elements can easily dominate over the elastic cross- 
section for a large range of energies. This was noted in ref . [9] and we expand upon 
this point here. 

Heavy elements, such as C", 0\ Si', Fe' , are not completely ionized but have 
their atomic inner shells filled. The total photoelectric cross-section (in units with 
h = c = 1) for the inner K shell mirror electrons of a mirror element with atomic 
number, Z, is given approximately by[28] 



16^2 
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me 
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(16) 



Evidently, the photoelectric cross-section decreases with mirror photon energy like 
{E'^)~'^/'^ and, of course, E'^ must be larger than the mirror electron binding energy 
of the particular element concerned. The contribution to the optical depth due to 
such inelastic scattering for 7' propagating out from the galactic center is 



Tis = / apEUA'dr 



~ J2 2poroapE 

A' 



rriA' 



(17) 



where ^a' is the proportion by mass of the mirror metal component. A' (e.g. A' — 
C ,0', Si', Fe', ...) and we have included a factor of two since there are two K shell 
mirror electrons. For illustrative purposes we have evaluated the total optical depth, 
including both elastic and inelastic scattering (the latter assumed dominated by K 
shell bound electron scattering as discussed above) for an example with a 2% metal 
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component with ^c' — Co' = ^Si' — ^Fe' — 0.005. The result is shown in figure 1. 
This figure indicates that for mirror photon energies 

0.4 keV keV (18) 

inelastic scattering of mirror photons can dominate over elastic scattering. Thus, 
it might actually be possible for supernovae to transfer a significant part of their 
energy to the halo in a fairly efficient manner (i.e. fsN ~ 0.03). 

Importantly, the optical depth is approximately independent of galaxy size since 
it scales as PqTq which is inferred to be constant [cf. Eq.(12)]. This means that the 
fraction of supernova energy transmitted to the halo is also roughly independent of 
galaxy size, which recall, was a necessary condition for the scaling relation, Eq.(ll), 
to hold. 
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Figure 1: Galactic optical depth, r, versus mirror photon energy (solid line). Also shown 
is the contribution to the optical depth due to elastic (Thomson) scattering, tes (dotted 
line). 



The scaling relation, Eq.(ll), was derived from only the very general condition 
that the total energy lost due to radiative cooling is replaced by the total energy 
input. Of course in any particular volume clement the energy input must match 
the energy output for a steady state configuration. It is useful to first consider a 
toy model, where we model the baryonic component as a point source whose energy 
output supports a spherical mirror dark matter halo. That is, we can assume a 
mirror photon luminosity, L, originating at r = 0. The energy going into a volume 
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element, dV — Airr^dr, assuming mirror radiation dominates the energy transport 
is 



dSin — FandV 



= -^andV (19) 

where a is the relevant mirror photon scattering cross-section [i.e. Thomson cross- 
section for E'^ ~ 30 keV and possibly a much larger photoelectric cross-section for 
E!^ in the range given in Eq.(18)]. We have assumed in this toy model that the 
optical depth, t = J andr ^ 1 so that the 7' flux, F, scales as ~ The energy 

going out of the same volume element is 

dE^t = KT)n''dV . (20) 

Matching d£in = d£out implies 

n^F^ = — — (21) 

A(T) 47rr2A(r) ' ^ ' 

The obtained n oc l/r^ behaviour also matches the condition from hydrostatic equi- 
librium of a pressure supported self gravitating spherical distribution with a common 
(i.e. independent of r) temperature, T[9]. Thus a self consistent 'toy model' emerges, 
except it is unphysical at r = 0. 

In the more realistic case there are two important perturbations to this picture. 
First the supernova energy sources are not all located at a point, r = 0, but a 
distribution extended over the disk scale length, r^j. This would be expected to 
soften the F oc l/r^ behaviour at small r ~ t-d. Secondly, as figure 1 suggests, 
t{E'^) 1 for a range of E'^. This effect can increase the proportion of supernova 
energy being absorbed in the inner region of the halo versus the outer region. If 
this effect dominates, then the inner (core) region will heat up relative to the outer 
halo region. If the temperature of the core does indeed rise then the mirror metal 
components can have their inner shells completely ionized. If this happens then these 
elments can much less effectively absorb supernova energy. There will ultimately 
be a balance at some temperature ~ Tc ~ i?^ ' . [Here E[^ 8 keV is the binding 
energy of the K shell mirror electrons for Fe'.] Thus we are lead to a picture where 
the inner region of the halo has temperature T ~ Tc and the outer region has a 
somewhat lower temperature given by Eq.(9). The outer region should be roughly 
isothermal and we expect pdm oc l/r^, while the rise in temperature of the inner 
region will result in a much flatter or even decreasing dark matter density as r — > 
(see discussion in the appendix). 

The transition radius between the lower temperature outer halo and hotter inner 
halo, which roughly defines the core radius Tq, will occur when (tjs) ~ 1- Here (r/5) 
is the appropriate averaging of tis{E'^) weighted by the 7' energy spectrum that 
supenovae ultimately produce. It might well be that the peak of this distribution is 
much greater than 10 keV (as indicated by the rough analogy with GRB's mentioned 
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earlier), in which case the faUing 7' spectrum in the relevant energy range E'^ ~ 30 
keV suggest that (tis) ~ tis{E'^/ ~ 10 keV). [Although the precise value of {tjs) is 
uncertain, the important point is that it is not expected to scale significantly with 
galaxy size.] Thus, these arguments suggest that the existence of the core in galaxies 
is related to the heating up of the inner region, which in this picture, should occur 
when Tjs{E' ~ 10 keV) ~ 1 or 



A' 



6 



ruA' 



apE{E'^ = 10 keV) 



\ 0.02 



(22) 



The line of reasoning culminating in the above equation provides a simple physical 
explanation for the observed scaling: Tq oc 1/po) indicated in Eq.(12). Moreover, in 
this picture the observed value of PqTq provides, in principle, a measurement of ^^e' 
which it turns out is compatible with direct detection experiments [3] and CMB/LSS 
constraints [8]. 

In conclusion, wc have considered galaxy structure within mirror dark matter - a 
dissipative and self-interacting dark matter candidate. For this type of dark matter, 
an energy source is needed to stabilize dark matter halos in spiral galaxies such 
as the Milky Way. Previously[9] it has been speculated that ordinary supernovae 
can supply the required energy if photon-mirror photon kinetic mixing of strength 
e ~ 10~^ exists. We have shown here that this argument motivates an approximate 
scaling relation, Rsn oc PqTq {Rsn is the supernova rate and poi ''"0 the dark matter 
central density and core radius). Interestingly, this scaling relation is consistent with 
the dark matter properties inferred from recent studies of spiral galaxies. We have 
also presented a line of reasoning, related to the details of supernova energy transport 
in the halo, which potentially explains why galaxies have a cored distribution and 
yields a second scaling relation, Tq oc 1/po- This relation is also in agreement 
with observations. Clearly more detailed studies are warranted to check the simple 
physical arguments presented here. 

Appendix - hydrostatic equilibrium 

Mirror dark matter, being self-interacting, forms a pressure supported halo. The 
pressure, P, can be determined from the hydrostatic equilibrium condition: 

dP 

^ - -P{r)9{9)^ 

= -mn(r)^ (23) 
r 

where p(r) = fhn{r), with rh being the mean mass of the mirror particles, e', H', He', .. 
in the plasma, n(r) total dark matter particle number density, g{r) is the local ac- 
celeration due to gravity. Here P — nT and Vrot{r) is the local rotational velocity. 
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For a spherically symmetric distribution, g (and hence v^^/r) can be related to the 
total mass density, ptotah via: 

g{r)^% f PtotaidV (24) 
r Jo 

As discussed in ref.[9], for an isothermal spherically symmetric halo, neglecting 
baryons, these equations have the solution T — ^fhv^g^ and p — X/r'^ = 2T^Gm) - 
Of course, the approximation of neglecting baryons could only be roughly true for 
the outer regions of the halo. 

Consider the case where the halo is not isothermal, but has a significantly hotter 
inner region. As briefly discussed already, this can be due to the efficient heating of 
the inner region due to the proximity of the supernova energy sources and significant 
opacity. The increasing temperature in the inner region means that the number 
density, n(r), should be much flatter than n{r) oc (or even decreasing as r ^ 0). 
If n(r) were constant in the inner region, then integrating Eq.(23) implies 

fro ^2 

P(0) - P(ro) ^rfin — dr (25) 
Jo r 

In the inner region of spiral galaxies, the rotational velocity, Vrotii^), is roughly 
constant until very near the galactic center where Vrot ~^ (see e.g.ref.[19]). It 
follows that 

T(0) - r(ro) ^ mvUro) (26) 

Given the isothermal result that T is related to the asymptotic Vrot value by: T = 
^fhv^g^, Eq.(26) indicates that a relatively modest increase in temperature in the 
inner region (less than a factor of two or so) would occur if n(r) were constant 
there. As discussed in the text, it seems reasonable to expect the temperature of 
the central region to significantly increase until the mirror metal components arc 
(almost) completely ionized, which suggests a central temperature around 1-10 keV. 
This is to be compared with the 'average' halo temperature, typically less than a 
keV, as indicated in Eq.(lO). We thus anticipate that the temperature difference 
between the central region and outer halo region will, typically, be somewhat larger 
than that suggested in Eq.(26). If hydrostatic equilibrium were to hold then this 
suggests that the number density, n{r), will likely be decreasing as r — )■ 0. Of 
course, this should not be a problem for observations (rotation curves), since the 
inner region of galaxies are well described by ordinary baryons [2 6]. 
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